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Abstract 

In this paper, we are concerned with nonparametric estimation of the multivariate 
regression function in the presence of right censored data. More precisely, we pro- 
pose a statistic that is shown to be asymptotically normally distributed under the 
additive assumption, and that could be used to test for additivity in the censored 
regression setting. 

Key words: censored regression; additive model; curse of dimensionality; 
nonparametric regression; marginal integration. 



1 Introduction and motivations 

A well known issue in nonparametric regression estimation is the so-called 
curse of dimensionality, i.e. the fact that the rate of convergence of non- 
parametric estimators dramatically decreases as the dimension of the co- 
variates increases (see, for instance, IStond (119821 )). To get round this issue, 
one common solution is to work under the additive assumption, i.e. the true 
regression function is assumed to be the sum of some lower dimension re- 
gression functions (typically, univariate or bivariate functions). But, this as- 
sump tion is strong and ha s therefore to be checked v i a one of the a v ailabl e 
tests (ICamlong-Viotl (l200lh . iGozalo and Lintonl (l200lh . ISperlich et all f liooih . 
Derbort et all (120021 )) before being used in practice. 



When the variable of interest is censored, several nonparametric estimators 

have been proposed for t he mu l tivariate regression f unction (see. e.g.jFan and Giibels 
( |l994J ). ICarbonez et all ( fl995h . IrCohler et a,l\ (120021 ) . iBrunel and Comtel (l200fih ). 
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By combini ng one of this 'initial' estimator with th e marginal integration 



method (see iNewevI (ll994T ). iLinton and Nielsen] ril995h ). estimates can be ob 



tamed under the additive assumption. In particular, iDebbarh and Viallon 
( 120071 ) made use of an initial Invers e Probability of Censor ing Weighted esti- 
mator (such as the one proposed by Carbonez et al. ( 19951 )). and established 
the uniform convergence rate for the corresponding additive estimator. How- 
ever, in this censored setting, no test for additivity has been proposed yet. 
That will be our concern here. Namely, we first exhibit a statistic evaluat- 
ing a weighted difference between the observations of the variable of interest 
and the estimator we derive via the marginal integration method. Then, this 
statistic is shown to be asymptotically normally distributed under the additive 
assumption. 



To build our estimators, and then our test statistic, the following notations 
are needed. Let (Y,C,X), (Y^C^Xi), (Y 2 , C 2 , X 2 ), ... be independent and 
identically distributed IxKx R -valued random variables. Here Y is the 
variable of interest, C the censoring variable and X = (X 1} X^) a vector 
of concomitant variables. In the right censorship model, the only available 
information on (Y, C) is given by {Z, S), with Z = min{Y, C} and 5 = 1 



Ie standing for the indicator function of the set E. As a matter of fact, 
observed sample is V n : = (X», Z if Si)i <i<n , for a given n > 1. 



{Y<C}, 

the 



Given a real measurable function tjj, our concern here is the regression function 
of ip(Y) evaluated at X = x, that is, 

m^x) = E (4>(Y) | X = x) , V x = (x 1; ...,x d ) G R d . (1.1) 

Under the traditional additive assumption, the regression function defined 
in (11.11) can be written as the sum of some (unknown) univariate regression 
functions mi, 

d 

J7ty(x) = m^ tadd (x) ■= n + mi(xi). (1.2) 

l=i 

In view of (II. 2p . the functions mi, as well as the constant term fi, are defined 
up to an additive constant. Therefore, we will work under the common iden- 
tifiability condition Emi{X{) = 0, for I = l,...,d. This condition implies that 



In the sequel, we set, for all t e R, F(t) = P{Y > t), G{t) = P(C > t) and 
H(t) = P{Z > t) the survival functions p ertaining to Y, C and Z r espectively 
Further denote by G n the Kaplan- Meier ( Kaplan and Meierl ( 1958 )) estimator 
of G, -.5cm 



G n {y)= n 

Ki<n 



NniZi 



N n {Zi) 



for all y > 0, with fa = I 



{Zi<y} 



1 - 5,. 



;i.s) 
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Here we defined N n 
were adopted. 



x 



\zi<x}, and the conventions Y[ = 1 and 0° = 1 



Consider the null hypothesis 
H : G Madd ■= {m : K 



L,m(x) = /i + ^m / (xi);E(m i (X;)) = 0}. 



Following the ideas of lHardle and Mammenl (119931 ) , ICamlong-Viotl (120011 ) and 
Gonzalez-Manteiga et al\ (120021 ). we denote by g some fixed weight function, 
by L a given kernel, i.e. a real measurable function integrating to 1, defined in 
R d and by (£ n ) n >i a sequence of positive real numbers. Further let m^ add be 
some estimator of under the additive assumption (II .2p (or, equivalently, 
under Ho). Now, let us consider the statistic 



rj-i-k 



,M -n i=l 



Siip(Zj) 



if>,add 



'X,; 



/n(X<) 



#(x)dx, 



;i.4) 



which is a natural estimator of the quantity E(E( S ^^ — m^ jCl(W (X))|X) 2 . 
Under a useful independence condition (see (CI) below), the latter quantity 
equals E(E(ip(Y) — m^, ia ^(X))|X) 2 , and then equals zero if and only if the 
hypothesis Hq is true. Moreover, in Theorem 12.11 below, this statistic is shown 
to be asymptotically normally distributed under Ho- Therefore, it could be 
useful to test for additivity in censored nonparametric regression. Properties 
of the corresponding test will be studied elsewhere. 



Now, we precise how m^ add may be constructed. Let K±, K 2 , K 3 and K, be 

kernels respectively defined in R, W 1 and W 1 . Further set f n the kernel 

estimator of /, with / denoting the density function of X. Namely, 



/»(x) 



n K j= i 



x, 



h, 



where (h n ) n >i is a given sequence of positive real numbers. Denote by (hj, n ) n >i, 
j = 1,2, two sequences of positive real numbers. To estimate the multivari- 
ate regression function defin ed in (11.111, the f ollowi n g Nadaraya- Watson type 
estimators can be used (see ICarbonez et all ( 119951 ) , iKohler et all (120021 ) and 
Jones et all fjl994j >). 

V^^)^ W ith W^(x) = ^fei, (1.5) 



and, for I = 1, 



■,d, 



m tp,n,l\ X ) 



yw l 



i=l 



^ )5 fw) with 



nh^nhz /n(Xi) 



;i.6) 
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where we set, for all x = (xi, ..,Xd) G R d and every I = 1, ...,d, x_z = (xi, .., 
%i-i,%i+i, •-, ^d)-To esti mate the addit ive c omponents, we use t he m arginal 
integration method (see Neweyl ( 1994 ) or Linton and Nielsen ( 1995 )). Let 
qi,...,qd be d given density functions. Then, setting g(x) = Y\f =1 qi(xi) and 
q-i(x.-i) = Uj^iqj(xj), we define 



r}i( x i) 



/ mw,(x)g_/(x_j)t2x_z - / m^(x)g(x)dx, Z = l,...,d, (1.7) 
in such a way that the two following equalities hold, 

r)i(xt) = mi(xi) - mi(z)qi(z)dz, l = l,...,d, (1.8) 

d „ 

m^(x) = V?7i(x/) + / mJz)q(z)dz. (1.9) 



In view of (11.81) and (11.91) . the functions r#, / = 1, turn out to be some 
additive components, and, from (11.61) and (11.71) . a natural estimator of the l-th 
component rji is given, for alH = 1, d, by 



[xi) = m£ { (x)g_ { (x_j)dx_,- / mjj, ,(x)g(x)dx. (1.10) 



From (ll.lOp . an estimator m^, add of the censored regression function can be 
deduced under the additive assumption (11.21) (or, equivalently, H ), 



d „ 
^V>,a^( X ) = S + / „ ^,n( X M x )dx. 



2 Hypotheses and Results 



These preliminaries being given, we introduce the assumptions to be made 
to state our results. First, consider the hypotheses pertaining to (Y, C, X). 
We suppose that (X, Y) has a joint density /x,y- Moreover, we impose the 
following conditions. 



(C.l) 
(C.2) 
(C.3) 
(C.4) 



C and (X, y) are independent. 
G is continuous. 

There exists a constant M < oo such that sup 0<t<T |^(t)| < M. 
m^, is a fc-times continuously differentiable function, k > 1, and 



sup 

X 



< oo; / 



c/. 



Remark 1 noteworthy that condition (C.l) is str onger than th e condi- 

tional independence of C and Y given X ; under which Weran U98f\) worked 
to build an estimator of the conditional survival function (see also \Dabrowska 



4 



( 199 A) ). Note, however, that the two assumptions coincide if C and X are 



independent. In other respect, to use Beran's local Kaplan-Meier estimator, 
the censoring has to be locally fair, that is P[C > t | X = x] > whenever 
P[Y > t | X = x] > 0. Here (see assumption (A)(m)(6) below), we essentially 
suppose that G(t) > whenever F(t) > 0, which is, on its turn, a weaker 
assumption. For a nice discussion on the difference between Beran's estimator 
and Inverse Pr o babili ty of Censoring Weighted type estimators, we refer to 
Carbonez et al. 



Denote by Ci, Cd, d compact intervals of M and set C = C\ x ... x C d . For 
every subset £ of R q , q > 1, and any a > 0, introduce the a-neighborhood 8 a 
of £, i.e. £ a = {x : ini ye£ \\x — y\\^q < a}, \\ • ||h« standing for the euclidian 
norm on M. q . 

We will work under the following regularity assumptions on / and fi, I — 
1, ...,d, fi denoting the density function of Xi. These functions are supposed 
to be continuous and we assume the existence of a constant a > such that 
the following assumptions hold, 

(F.l) : Vxi G Cf, fi(xi) > 0, 1 = 1, d, and Vx G C a /(x) > 0. 
(F.2) : / is fc'-times continuously differentiable on C a ,k' > kd. 

Regarding the weight function g, we will assume that the condition (G.l) 
below is satisfied. 

(G.l) : g is an indicator function with compact support included in C. 



Remark 2.1 Assumption (G.l) is made here to avoid technical issues in the 
derivation of our results. Moreover, it is not restrictive since g is a given 
weight function. That being sai d, this assumption cou ld be relaxed (see, e.g., 



Gonzdlez-Manteiga et al.l $2002) . \ Camlong-Vioi $20 Oft) ). 



The kernels L, K, and K% defined in IR d , K\ defined in R and K<i defined in 
IR^ 1 , are assumed to be continuous, compactly supported and integrating to 
1. Moreover, we suppose that, 

(K.l) : K\ is Lipschitz; 

[K.I) '■ K\ and K 3 are of order k, and K is of order k' . 

In addition, we impose the following assumptions on the integrating density 
functions q_i and qi, I = 1, d. 

(Q.l) : q-i is bounded and continuous, / = 1, ...,d. 

(Q.2) : qi has k + 1 continuous and bounded derivatives, I = 1, d. 

Turning our attention to the smoothing parameters £ n , h n and hj >n , j = 1,2, 
we will work under the conditions below. 
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'lograN i/(2fc'+d) 
(id.l) : /i n = ci ( ) , for a given < C\ < oo. 



'logn\ V(2fe+i) 

(i?.2) : /ii in = C2 ( I , for a given < c 2 < oo and /i2,n = o(l). 



(F.3) : n{\ogn/n) k/{2k+1) t d J 2 -> and <! 



oo. 



As mentioned in iGross and Lail (119961 ) . functionals of the (conditional) law 
can generally not be estimated on the complete support when the variable 
of interest is right-censored. Accordingly, we will work under the assumption 
(A) that will be said to hold if either (A)(z) or (A)(ii) below holds. Denote 
by Tl = sup{t : L(t) > 0} the upper endpoint of the distribution of a random 
variable with right continuous survival function L. 

(A) (2) There exists a r < Th such that if) = on (r , oo). 

(A)(m) (a) For a given k/(2k + l)<p< 1/2, fin F-P/^dG\ < oo; 
(6) T F <T G ; 

(c) n 2p ~ 1 h ln \ log(/ty n )| —> oo, as n — > oo, for every Z = 1, d. 

It is noteworthy that assumption (A)(m) allows for considering the estimation 
of the "classical" regression function, which corresponds to the choice ip(y) = 
y. On the other hand, normality for estimators of functionals such as the 
conditional distribution function P(Y < 7o|X) can be obtained under weaker 
conditions, when restricting ourselves to t < Th- 

To state our result, some additional notations are needed. Set e» = ^j$r§\ — 
m^,(Xj) andag(x) = i£(e 2 |Xj = x). Further introduce B = [J o"q(u)/ _1 (u)(7(u) 
du] x / L 2 (t)dt and V = 2 [/(a-g(u) 2 /~ 2 (u)0 2 (u)du] x / [/ L(t)L(t - r)dt] 2 dr. 



Theorem 2.1 Assume the conditions (A), ((7.1-4), (F.l-2), ((7.1), (K.l-2), 
(Q-l-2) and (#.1-3) /io/d. Then, under the null hypothesis H , we have, 



V 



Af(0, 1), as n 



oo. 



3 Proof of Theorem [2~T1 



Here, we present the proof of Theorem 12. II in the case where (A)(z) holds. The 
case where ( A)(ii) holds follows from similar ar guments (especially replacing 
the result of iFoldes and Reitol (119811 ) by that of iGu and Lail (Il990l ) or that of 
Chen and Lol ril997t )): details are then omitted. 
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We will make frequent use of the following lemma, which was established in 
Debbarh and Viallonl d2007h . 



Lemma 3.1 Assume H holds. Then, under the conditions (A), ((7.1-4), 
(F.l-2), (if. 1-2), (Q.l-2) and (H.l-2), we have, with probability one, 



sup dd {x) -m^xj 

xGC 



o 



logn 



n 



2fc+l 



(3.1) 



We wi ll also make frequent use of the following result, due to lFoldes and Rejto 
(jl98lh . 



For all r' < T H , sup \G* n (y) - G(y)\ = (7((loglogn/n) 1/2 ) =: p n (3.2) 

y<r' 

From this last result, we especially get the following type of approximations. 
Set e* = gg} -m^(Xi). Then, from (A)(i), ((7.2), ((7.3) and ([372]), we have, 
almost surely as n — > oo, 

e^ = ei + (7(p n ). (3.3) 
Now, recalling the definition (11.41) of T*, we have 



1 " ,x-XA/m v ,(X J )-mV (M (X i ) + eJ 



-"*ti i=l n 

Consider the quantity 



g(x)dx. 



1 " /x - ^m^(Xi) - rn^ >add (Xi) + e\ 

"^n i=l 



*\n2 



4 / V /(x,) 

with fsee ICamlong-Viotl (l200lJ )). 



g(x)dx. 



T i* 

I n,l 



J n,2 



J n,3 



J n,4 



J n,5 



Rd n 2£2d * 



1 A 

1 E^f^W^^ .g^Lcfac, 



n //(X i )/(X j 



<fe V 4 A /(x,) 



<J i=l 



flf(x)dx, 



0* 



>^(Xi) -m^ add (Xj))e*-^— — dx 



/x - XA /x - XA fm^Xi) - m^^X*, 
2>( — ) L \ — ) ( TrsFTJrsr \ )e,^(x)dx. 



t n ) V /(X i )/(X i ) 
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By (F.l- 2), (H.l) and (if. 2) , it ho lds that, almost surely as n — > oo (see for 



instance 



sup 

xec 



Aneo-Nze and Riosl (l2000h ). 



/n(x) - /(x) 



O 



'log ?2 



hi P 



(3.4) 



in such a way that T* = T^*{l+0[(n 1 h n d logn) 1 ^ 2 ]} almost surely as n — ► oo. 
Therefore, to achieve the proof of Theorem I2.1[ it is sufficient to establish 
(Q-O below. 



n£ d X:=B + o p (£ d n /2 l 



n- 



(d/2 



( T n*2 ~ ET n,2) 
V 

?£, = o p {n-% d '\ 
< 5 = o p {n-X d ' 2 )- 



N{0,1), 



(3.5) 
(3.6) 

(3.7) 
(3.8) 
(3.9) 



Proof of ( 13. 5 j) : Set o"q 2 (x) = E(e* 2 |Xj = x). Using a conditioning argument, 
it is straightforward that 



ET nA 



n£l d 



gg(v) 
/(v) 



x — V 



dvg(pc)dx. 



Moreover, since g is an indicator function with compact support included in 
C, we obtain, for n large enough, that 

ET n,i = ^a j E(ef |Xx = x)^(x)r 1 (x)dx / L 2 (u)du + oQ, 
and then, wo arguments similar to those used to derive fl3.3|) . 



ET n,i = —M B + c '(pn), asn^oo. 



Turning our attention to the variance of T^* l5 we can write, 



VarT^ JLvar(I), where I = J e\ 



2 L 2 



' x ~ x j\ #( x ) 



/ 2 (X 



■dx. (3.10) 



But, using once again the arguments used to show (13. 3p . along with the facts 
that G{t) > 0, if) and g are bounded and L is compactly supported, it holds 
that, as n — > oo, 



£(i 2 



S(ef|Xi=v) /L 



x-v\ g{x 



/ 2 (v 



-cfx /(v)dv, 



0(O(l + ^(Pn)). 



(3.11) 
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By flXTUl) and (EOT]) , it follows that Var T^ = 0(n- 3 £- 2d ). Using the Bienayme- 
Tchebytchev inequality, we infer that, for all e > 0, P[n£ d [ 2 \T^ 1 — ET^\\ > 
e) < 6~ 2 nX x VarT^ = 0((n£)-i). Thus, = (ni^B + 'o p [(n£ d /)- 1 }, 
which is (13.51) . 



Proof of (SSD: For 1 < i < n, set Q = (X,, e*) and = L[(x - Xi)/4]- 
Further introduce 2 = ^ 2d fX)i<j M «( x )' u i( x ) e i e i5 , ( x )/ _1 (Xj)/ _1 (Xj)c/x. Note 
that, in view of (13. 2p and (13.41) . the dominated convergence theorem ensures 
that T}* 2 — T* 2 + 0(p n ) almost surely as n — > oo, with T* 2 := 2n~ 2 T^ 2 . 
To est ablish 03.61) . we will make use of a central limit theorem for U-statistics 
duetoid3([li9). Set M n (Ci,0) = C d / u i (x) Mj (x))e i e^(x))/- 1 (X l )/- 1 (X j )rfx 
and iV n (u, v) = E(M n (Ci, u)M n (d, v)). To apply Hall's theorem to T^ 2 , the 
conditions [Tl], [T2] and [T3] below must be verified. 



[21] £?{M n (Ci,C2)|Ci} = 0. 
[T 2 ] E{M2(d,C 2 )}<oo. 

[Ik] \E{N 2 n ((u(2)} + n-'E{M*(( l ,(2)}\ / \E{M 2 n {Ci, C 2 )}| 2 - 0. 

[Ti] is readily satisfied by making use of conditioning arguments. Moreover, 
arguing as before, the statement [T 2 ] follows from routine analysis. To establish 
[T 3 ], it is sufficient to prove the results (I3.12p . (13.131) and (13.141) below. 



E(N 2 (( 1 ,(2)) = 0(l) asn^oo, 
E(M 2 (( 1 ,(2))=£^ + o(t n 

£(M n 4 (Ci,C 2 ))=0" 



as n 



oo, 



-3d 



as n 



oo. 



(3.12) 
(3.13) 
(3.14) 



Proof of \3.12i) : Denote by /x, e the joint density of (X, e) (the existence of 
which being ensured by the assumption (C.l), since (X, Y) is supposed to 
have a joint density). It holds that 



E 



xL 



N 2 n ((i,C2 



L[ Xi : Vi WyU /X2 ' V| 



x 2 



Hi) 11 



XlJ #(x 2 J 



in ) / 2 ( Vl ) / 2 (V 2 ) 



/x,e(vi, ei)ele 2 e 3 dvideid]Ci(hc2 



de2dezdrjidrj2. 



Using classical changes of variables, together with the assumption (C.l) and 
the fact that L is compactly supported, (13.121) is straightforward. 

Proof of (E23). Set a 2 = JJ M 2 (oo u ^f^Mfr^du^,. Then, 



1 ^/ Xi - Vi \ / Xi - v 2 

£i d V 



L( x 2 - Vl x /X2 -v 2 



"'• ~ 11 1 '' ' ' ' 1 -dxidx.2(tvidv2deide2. 



xe 2 e 2 / Xl , ei (vi,ei)/x 2 ,e 2 (v2,e2) 



/ 2 (v 2 ) / 2 K) 
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Next, noting that / e 2 /x,e(vi, e\)de\ = £?(e||Xj = vi)/(vi), it follows that 



o: 



x = v 1 )/(v 1 )i?(e 4 2 |X J = v 2 )/(v 2 ) 9 f ^ 2 \ ^. d^d^d^dv,. 

r(v 2 ) r(vi) 

Using the changes of variables, yi = (xi —Vi)/£ n , Y2 = (x 2 — vi)/£ n and rj = 
(vi — v 2 )/4i, along with the continuity of / and the dominated convergence 
theorem, we get, 



J ... J ^L( yi )L(y 2 )L( yi - n)L(y a - r 1 )E(e l 2 |X, = Vl ) 

v!?/ 2i Y . o + yi£ n )g(vi + yi4) , , , , 

x£(e, X 3 - = vi + rx4) — — — ■ —dy 1 dy 2 dr l dv 1 , 

/(vi) /(vi + ri4) 



= C*/[/ L(t)L(t - r)dt] rfr JiE^Xi = r)) 2 |^dr + o(C 
which, recalling the definition of V, is (13.13!) . 



Proof of \3.1\ ): Arguing as before (see also ICamlong-Viot we can 

show that, for a given C < oo, 



|£(M n 4 (0,0))l <S 



L( x 3 -v 1 X /X 4 - Vl X / Xl -V 2 X /X2 - v 2 



L( -X3-V 2 y/X 4 -V 2 



g(x 3 ) g(x 4 ) 
/ 2 (vi) / 2 (v 2 

o(C d ). 



/(vi) /(v 2 ) 



■ dv i (i v 2 <ixi dx 2 c?X3 cfocj 



Combining (13.121) . (13.131) and ( 13. 14ft . it is readily shown that [T3] holds. Then, 
Hall's Theorem can be applied to T^ 2 . Namely, since EM n (Q, Q) = 0, we have 
^/2{na n )^ l T^ 2 — > jV(0, 1). Recalling that T* 2 = n 2 T^' 2 /2, we deduce, from 
(I3TT3D . that {p£% 2 )T^/y/V -> JV(0,1). Slutsky's Theorem is now sufficient 
to conclude to (13.61) . since, as already mentioned, T^ 2 = T^ 2 + 0(p n ) almost 
surely as n — > oo . 



Proof of (15771) : By (G.l), 

T n,3 = SU P l«^( X ) - ™^( X ) 



x- X,- 



g(x)dx a.s. 
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Since, under H , = m^^dd G -Madd, we can apply the result of Lemma 
13.11 This latter, when combined with, successively, the boundedness of /, 
the dominated convergence theorem, Bochner's theorem and the fact that g 
is compactly supported, yields T** 3 = 0f(logn/n) 2fc / 2fc+1 ) almost surely as 

n — > oo. Thus, under the assumption (H.3), T^* 3 = o p (n~ 1 £~ d ^ 2 ) as n — ► oo. 
Proof of ( 13.81) : First consider the mean of T^ 4 . By (G.l), it holds that 



I ET^ A | < sup | (x) - m^ add (x) | ^ ^ 
Xg(x)dx.du. 



E[ e*\Xi = u)i/ 



x — u 



Next, under the assumptions (if. 1-2-3), using, successively, the assumption 
(G.l), the dominated convergence theorem, the equality (13.31) . Bochner's theo- 
rem and Lemma [37TJ it can be shown that |_ET^ 4 | = O ((ni^) -1 (log n/n) k '^ 2k+1 ^) 
= oipr 1 ^! 2 ). Turning our attention to the variance of T^* 4 , and arguing as 
before, we get 



Var Tl a = - „, , Var 

TV 



X - X A/._. / V N ^ / V ^ #( x ) 



e\L 2 — — i {m+VQ - m* (X^i^dx 



/ 2 (X.) 



/ 1 / logn \ 
U 3 ^A n / 



W n \ 2A/(2A+1) 



An application of Bienayme-Tchebychev's inequality leads to T* 4 — ET\ A 
Opfa- 1 ^*' 2 ), which implies (I3T51) . since |£T** 4 | = ofa -1 ^/ 2 ). 

Proof of (13.91) : Arguing as before, we infer that, ultimately asn^ oo, 
//logn\ fc /( 2fc+1 )\ / 1 /logn\ 2fc /( 2fc+1 ) 

Therefore, Bienayme-Tchebychev's inequality leads to ( 13. 91) . 
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